TeVeS gets caught on caustics 
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TeVeS uses a dynamical vector field with timelike unit norm constraint to specify a preferred local 
frame. When matter moves slowly in this frame - the so-called quasi-static regime - Modified Newto- 
nian Dynamics (MoND) results. Theories with such vectors (such as Einstein-Aether) are prone to 
the vector dynamics forming singularities which render their classical evolution problematic. Here 
we analyse the dynamics of the vector in TeVeS in various situations. We begin by analytically show- 
ing that the vacuum solution of TeVeS forms caustic singularities under a large class of physically 
reasonably initial perturbations. This shows the classical evolution of TeVeS appears problematic 
in the absence of matter. We then consider matter by investigating black hole solutions. We find 
large classes of new black hole solutions with static geometries where the curves generated by the 
vector field are attracted to the black hole and may form caustics. We go on to consider the full 
dynamics with matter by numerically simulating, assuming spherical symmetry, the gravitational 
collapse of a scalar, and the evolution of an initially nearly static boson star. We find that in both 
cases our initial data evolves so that the vector field develops caustic singularities on a time scale of 
order the gravitational in-fall time. Having shown singularity formation is generic with or without 
matter, Bekenstein's original formulation of TeVeS appears dynamically problematic. We argue that 
by modifying the vector field kinetic terms to the more general form used by Einstein- Aether this 
problem may be avoided. 

PACS numbers: 04.50.Kd,04.70.-s,04.25.D-,04.20.DW,04.70.Bw,97.60.Lf 



I. INTRODUCTION 

For many years the existence of dark matter has been 
postulated to reconcile a number of astrophysical and 
cosmological observations with our understanding of the 
laws of gravitation. Dark matter was originally intro- 
duced decades ago to explain the discrepancy between 
the rotation velocities of stars in the outskirts of galaxies 
and that predicted by the mass inferred from the amount 
of visible mass in the galaxy. 

The success of the dark matter paradigm extends be- 
yond galactic scales to cluster and indeed cosmological 
scales. Today we know that best fit models of struc- 
ture formation apparently require a dark matter fraction 
much larger than the known baryon content of the uni- 
verse to drive the growth of structure from kpc through to 
Gpc scales. The potential wells provided by a cold, dark 
matter (CDM) component also reconcile the amplitude 
of the acoustic peaks observed in the Cosmic Microwave 
Background (CMB) angular power spectrum [l|, [3, d, 3 
with the known baryon content. 

Dark matter also provides a simple explanation for the 
observed lensing of background galaxies by clusters along 
the line of sight. Recently the combination of optical, x- 
ray and lensing observations of the bullet cluster have 
yielded the most direct evidence to date in support of 
the picture where the gravitational mass of clusters is 
dominated by a dark matter component 



Taken as a whole the growing wealth of observa- 
tions points clearly to a concordance ACDM cosmologi- 
cal model with a significant fraction of the critical energy 
density made up of CDM. The dark matter paradigm has 
stood the test of time remarkably well but significant 
questions remain. Many candidates for a dark matter 
particle exist ranging from massive neutrinos to more ex- 
otic weakly interacting extensions to the standard model. 
However dark matter has yet to be detected directly in 
the laboratory or indirectly possibly through the 7-ray 
signature of its decay in the centre of galaxies (this is 
required to avoid the concentration of dark matter ob- 
served in numerical simulations). 

For these reasons an alternative approach to adding 
a dark matter component has been to consider whether 
the discrepancies between observations and general rel- 
ativity in the low acceleration regime are an indication 
of the failure of the_theory itself. This was the approach 
taken by Milgrom [6] who proposed a phenomenological 
modification to the acceleration equation which seems to 
fit well galactic rotation curves without the addition of 
any dark matter: 



(|a|/ao)a= -V«>, 



(1) 
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where $ is the Newtonian potential, ^{x) is an arbitrary 
function with limits such that ^{x) ^ 1 in the strong 
acceleration regime {x ^ 1). The constant ~ 10^^° 
m s~^ determines the acceleration scale below which the 
Modified Newtonian Dynamics or MoND becomes rele- 
vant, and the above acceleration law receives non-linear 
corrections. While is potentially a free (monotonic) 
function, only the limits where its argument goes to zero 
or infinity affect the astrophysical phenomenology. 
MoND has been successful in fitting the anomalous 
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accelerations observed in galaxies and clusters (see e.g. 
01 for a recent review). It also successfully predicts the 
TuUy-Fisher relation correlating the luminosity of galax- 
ies to the fourth power of the rotation velocity. However 
it remains a phenomenological modification of gravity 
with no underlying relativistic theory. In addition the 
simplest theory based on a modified action that repro- 
duces Eq[l] depends explicitly on coordinates and breaks 
conservation laws. The fact that MoND has no underly- 
ing covariant theory has restricted its application for the 
purpose of comparison with other astrophysical and cos- 
mological observations. For example, lensing predictions 
in both the strong and weak regimes cannot be formu- 
lated in MoND and this has left it unable to answer the 
criticism stemming from lensing mass reconstructions of 
galactic and cluster profiles which seem to suggest the 
existence of dark matter halos. 

Recently however, Bekenstein Q, has put forward 
TeVeS, a relativistic theory of gravity which reduces to 
MoND in the weak acceleration limit. In TeVeS the mat- 
ter sector lives on a matter frame (MF) metric which 
maps 'disformally' to a second, Einstein or gravitational, 
frame (EF) metric via a dynamical scalar field (p and 
a dynamical vector field A. The addition of a scalar 
and vector degree of freedom are behind the name 'Ten- 
sor, Vector, and Scalar', or TeVeS theory. TeVeS builds 
on previous attempts to obtain a relativistic version of 
MoND which suffered from a number of inconsistencies 
involving the acausal propagation of physical degrees of 
freedom ^, . TeVeS however was shown to be a fully 
causal theory for positive values of the additional scalar 
field. 

The original motivation behind TeVeS was to build a 
theory with a fully consistent action which recovers the 
MoND behaviour in the weak acceleration limit. How- 
ever, given it is a relativistic, metric theory of gravity 
and matter, it can do much more. In TeVeS it is possible 
to calculate geodesies in the presence of a matter sources 
which leads to lensing predictions [ll| . It is also possible 
to show that it is compatible with the basic background 
cosmological observations such as age and distance mea- 
sure observations Q. The full framework of relativistic 
perturbation theory can be developed in TeVeS which 
makes comparison to the perturbed universe possible. 
Already the first calculations in this area have shown 
that the theory may be reconciled with CMB and Large 
Scale Structure (LSS) observations [H, d [H, [Hi , albeit 
with some fine tuning of the model ingredients. Attempts 
have also been made to explain the bullet cluster results 
within the TeVeS framework p^ . 

For TeVeS to be a successful theory it must also be 
shown to be consistent, and agree with observations, in 
the strong gravity regime. In exploring this end of the 
theory the potential is that it could be compared to as- 
trophysical observations of compact objects such as neu- 
tron stars and black holes or at the solar system level 
with post Newtonian (ppn) corrections to planetary or- 
bits 



In order to have a modification of gravity dependent 
on acceleration, one must have a reference frame in which 
to measure that acceleration. The vector in TeVeS dy- 
namically selects that reference frame, spontaneously 
breaking Lorentz invariance since it is constrained to 
have unit timelike norm. All types of matter see the 
same distorted metric so adding a preferred frame is not 
in conflict with weak equivalence principle tests. Only 
tests of gravitational dynamics can constrain the theory. 
MoND is recovered from TeVeS when matter moves non- 
relativistically in the frame defined by the TeVeS vector, 
which has been termed the 'quasi-static' regime. The 
purpose of our work is to argue that this quasi-static 
regime will typically only exist for a short period of time, 
of order the gravitational in-fall time, after which the 
vector field develops a singularity and the theory can- 
not be classically evolved any further. Hence TeVeS even 
classically is dynamically sick in practice and recovery of 
MoND or even GR is impossible. Indeed here for simplic- 
ity we will focus on the large acceleration regime relevant 
on small scales (e.g. within the solar system) where the 
in-fall time scales are shortest. 

This singular vector field behaviour is analogous to 
that in other modified gravity theories such as Einstcin- 
jEthcr theory [19, 20,] and Ghost Condensation .21, 2^]. 
Einstein- iEther theory is much simpler than TeVeS, be- 
ing simply Einstein gravity modified by adding a vector 
field, again with timelike unit norm constraint. The vec- 
tor action is taken to be more general than that in TeVeS 
where it is simply that of a Maxwell field, but one may 
choose them to be the same. In this case (actually a 
theory written down earlier [2^), it is easy to show that 
the vector field generically develops singularities; classes 
of solutions exist where the integral curves of the vector 
are timelike geodesies moving in the spacetime geometry 
created by the matter. These geodesies fall into gravity 
potential wells and meet, and when they do so, the flow 
they define develops caustic singularities [l^ . The vector 
field at these points becomes singular. It is for this rea- 
son that the Einstein- jEther literature focuses on other 
choices of the vector action than Maxwell type. Indeed 
while the ghost condensation theory has no vector, it is 
the integral curves of the gradient of the ghost scalar that 
form caustics. Since TeVeS is a considerably more com- 
plicated theory than Einstein- iEther, with complicated 
coupling of its vector and scalar to the matter, the vector 
behaviour and in particular whether it forms singularities 
could be very different. Our key result is that while in 
detail the dynamics of the vector is clearly different, it is 
still subject to the same singularity development as the 
Maxwell case of Einstein- iEther. However, we can play 
the same game as in Einstein- iEther theory, and by tak- 
ing more general vector kinetic terms, we may avoid this 
behaviour, and as we show later, we still recover MoND 
for quasi-static systems. 

The paper is organised as follows. In section [L^ we re- 
view TeVeS theory and the field equations derived from 
the TeVeS action. We review the relation between TeVeS 
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and the Einstein- iEther theory in section II Bl where we 
introduce the problem of vector field caustic singularity 
formation. In section|ll]we begin by showing analytically 
that in the absence of matter the TeVeS vacuum forms 
singularities under evolution of a class of initial perturba- 
tions. These are physically reasonable perturbations and 
we are able to precisely characterize the condition for 
the initial data developing to a singularity. We find that 
the condition is generic within the class of perturbations. 
Already this analysis indicates a fundamental problem 
with the TeVeS dynamics. However, phcnomcnologically 
it is the inclusion of matter in TeVeS that is of key inter- 
est, and hence we proceed to study whether such inclu- 
sion ameliorates or worsens the problem. We begin our 
study of matter in TeVeS in section IIIII with an analytic 
treatment of black holes. Previous static black hole solu- 
tions of Giannios [l3| and Sagi and Bekenstein have 
a static vector field aligned with Killing time. However 
we find a large class of new solutions where the geom- 
etry is static but the vector field is in general dynamic. 
There exists a family of stationary solutions where the 
vector falls into the future horizon, but also dynamic so- 
lutions where caustics may form in the exterior of the 
black hole. The black hole provides a focus for the curves 
of the vector field and hence indicates that singularity 
formation is likely to be enhanced by compact matter 
sources. We then proceed in section ITVl with our study of 
matter by performing full numerical simulation of scalar 
collapse and the evolution of an initially quasi-static bo- 
son star. In both cases the evolution ends at a caustic, 
in the former outside of an apparent horizon, and in the 
latter near the surface of the star. Having given evidence 
that the dynamics of TeVeS is too pathological to provide 
a relativistic setting for the phenomenological theory of 
MoND, in section |V] we outline a modified TeVeS the- 
ory which may not suffer from caustic formation in the 
vector field. We explicitly show that MoND is again re- 
produced in the appropriate Newtonian limit. We argue 
that this modification is likely to affect many detailed 
phenomenological studies of TeVeS should therefore be 
included. We conclude with a summary and brief discus- 
sion of our main results in section IVTl 



A. TeVeS Action and Field Equations 

TeVeS is constructed using two metrics, the matter 
frame (MF) metric g, and the Einstein frame (EF) met- 
ric g. The two metrics are related through 'disformal 
relations' involving the extra scalar and vector fields 
and A„. 
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- 2A„A0sinh2(?!), (2) 
2A"A'^sinh20. (3) 



The total action S governing the dynamics in TeVeS can 
be split into separate components S — Sg^y -1-5*8-1- Sm, 
where 
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){~gflH^x, 



(4) 



where g is the determinant of the EF metric, i? is the 
scalar curvature, G is the gravitational constant and 
F^v = \tL,v\ = A^^v - A„^f^. The Lagrange multiplier 
A enforces the timelike, unit norm constraint on the vec- 
tor field 



g'^'^A^A, = -1. 
The scalar field action is given by 



(5) 
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{-gY'^d^x, (6) 



where ct is a non-dynamical scalar field and J-{KGa^) is 
a dimensionless function whose behaviour is determined 
by requiring GR and MoND to be recovered in the ap- 
propriate dynamical limits. TcVcS introduces three new 
parameters; the two dimensionless constants k and K 
and a third parameter / with units of length. 

Finally the matter action S„i is built using the MF 
metric as 



Srn= / i[5,X^,9x1 {--gy^'d^ 



(7) 



for a collection of matter fields x^- Thus all matter fields 
are coupled to the same MF metric and test particles fol- 
low the same geodesies. This ensures the weak equiva- 
lence principle is satisfied and the theory is not in confiict 
with fifth force measurements. 

Varying with respect to g, and recalling g = g (g. A, (j)) 
gives the Einstein equations 

G„;3 =87rG[f'„^ + (l-e-4*)A^f^(„A^)+T„/3]+e„0, (8) 

where i),(a^,/3) = Tf^aAp +ff_,pAa, Gap is the Einstein 
tensor, Taj3 is the energy momentum tensor of the matter 
components defined in terms of the MF metric cjap and 
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XA^Ap. (10) 



Variation with respect to the scalar field a yields a rela- 
tion between a and (jj^a involving J-'. The specific choice 
of J- determines the exact behaviour of the theory in the 
weak acceleration regime and is relevant for the MoND 
and cosmological behaviour of TeVeS. The regime of in- 
terest for this work is one where the acceleration is much 
stronger than the MoND acceleration scale qq. In this 
case the MoND function /i(|a|/ao) 1, which is equiva- 
lent to a limit on the argument of the free function J-', 
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(11) 



For any suitable function choice, diverges logarithmi- 
cally in (/i — 1) in this limit. The contribution of J- to 
the field equations ([5]) however is suppressed by a factor 
(/i — 1) relative to other terms, and so when /x ~ 1 it 
may be neglected [1, US]- Thus our results will be 
insensitive to any particular choice of J- and we drop the 
term in the following. Finally variation with respect to 
the scalar gives 



[ig'^^-A-A^)^,^],^^^ 
kG [g'^f^ + {l + e-^'^)A''A^]f^p, 
and for the vector we have 



(12) 
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A^<j>.,pg^''cl), 



87rG(l - e~'^''')g°'^'A'^f^p. 



(13) 



As stated here TeVeS is a classical phenomenological 
theory. The somewhat Baroque form for the Lagrangian 
leads to the obvious concern that the theory is not sta- 
ble to quantum corrections. Attempts have been made 
to study a UV origin from String theory [H, |2^, and 
there are certainly many interesting questions in these 
directions which we do not consider here. 

It is also worth mentioning that the TeVeS theory itself 
been generalized by various authors [3, 113, HI, [H, H^l 
and it would be interesting to consider the formation of 
caustics which we study here in these modified versions 
of the theory. 



B. ^ther theory, its relation to TeVeS, and 
problems with its vector field dynamics 



Another theory of aether field dynamics is Einstein- 



jEther theory 
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an effective field theory designed to 



investigate the effects of Lorentz violation in a fully co- 
variant setting. It has the action. 
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R + K'^^.Vo.A^^VpA- 



\{A^ + l)\{-gf/'^ + / CnatteAg]. 



(14) 



where K"^^^ provides the most general kinetic term for 
A which is diffeomorphism invariant, quadratic in deriva- 
tives and (preemptively) consistent with the A^ = —1 
constraint. Specifically, 
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a/3 



cig'^'^gt^u- 
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'C4 



A-A^g,,. (15) 



This kinetic term is the usual Maxwell case when c+ = 
ci + C3 = 0, c4 = c2 = 0, c_ = Ci — C3 < 0. Einstein- 
iEther theory is actually a truncation of TeVeS in the 
absence of matter, where we may consistently set the 
scalar to zero and then c_ = ~2K - however obviously 
phenomenologically this is not the regime of interest for 
TeVeS where the coupling to matter and the non-zero 
scalar are crucial. 

Following Jacobson & Mattingly [l^ it is easy to see 
that the Maxwell case of Einstein- jEther is pathological. 
To any solution of Einstein gravity coupled to matter, we 
may simply add a vector field obeying the equations. 



Fa. = 0, 



A^ = -l, 



(16) 



and this will then solve the full Einstein- iEther equations 
for that matter since the vector and constraint contribute 
nothing to the stress energy. Note that the vacuum, the 
Minkowski geometry and A^^ = (dt)^, is in this class of 
solutions. Generally the solution is given by. 



A,, 



{dxf 



-1, 



(17) 



where the latter equation is a p.d.e., first order in time. 



dtX = 



1 



(-3") 



(18) 

with i = 1, . . . , 3, which can evolve x in time from an 
initial Cauchy surface. We have taken the choice of root 
above since we wish A'^ to be a future directed vector 
field. Hence the data for the solutions can be charac- 
terized by the function x(i = 0,x). Now, A''^/ ^A^^ = 
A^V^Ajy = using both relations in (fT6|l . Hence integral 
curves of the vector field A are simply timelike geodesies. 

Suppose we consider a static star as a matter source. 
Then the solution above will have families of vector fields 
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with different initial directions, but all will have inte- 
gral curves that fall in towards the gravitational poten- 
tial well and will meet each other in a timescale of order 
the gravitational in-fall time. Since these are integral 
curves of the vector field, when they meet they result 
in a caustic singularity where the value of the vector is 
ill-defined. Indeed, even in the absence of matter and 
with the Minkowski spacetime geometry it is possible to 
have singular behaviour. Such a solution is illustrated in 
figure [TJ 

This simple argument shows for solutions with F^^, — 
that caustic singularities generically occur in the presence 
of gravitational potential wells. It seems reasonable that 
singularities will also occur in solutions where Ffj,i, ^ 0. 
Whilst there are no general arguments, in specific cases 
with Fu^y ^ singularity formation has been shown by 
Clayton [2(| |. 

It is for these reasons that the Einstein-jEthcr litera- 
ture does not consider the Maxwell vector kinetic term 
[3l| . Interestingly there is little rigorous understanding 
for what choices of parameters Ci do give well behaved 
vector dynamics. Clearly a caustic singularity is signalled 
by the divergence of the vector field becoming infinite. 
Hence it is expected that by adding the term C2 appro- 
priately, which directly energetically weights this diver- 
gence, one can dynamically suppress singularities. For 
example, taking the Maxwell case together with the ad- 
ditional term C2, we obtain a vector equation of motion, 

[81 + A^A^) {d.F''^ -t- C2d^{d ■ A)\ = 0. (19) 

We note that when C2 = 1 the equation of motion is essen- 
tially the wave equation, and hence in a regular geometry 
we would certainly not expect singular behaviour. How 
large C2 should be to avoid singularities is an interesting 
open problem. 

Ignoring gravity and matter, we simply plot a vector 
solution to the above equation in 1-1-1 flat space in fig- 
urefl] and [H The left frame is for C2 = 0, the right 
for C2 = 0.2, and both have the same initial data which 
satisfies = 0. Note that the right frame cannot be 
lifted simply to a solution of Einstein- iEther since for 
C2 = 0.2, F^j_y will not remain zero, and the vector neces- 
sarily contributes to the gravitational stress tensor. How- 
ever, clearly by eye we see a change in behaviour, where 
the small amount of positive C2 avoids the caustic, lead- 
ing to an asymptotic vector solution which is aligned with 
time. 



II. VECTOR FIELD DYNAMICS IN THE 
ABSENCE OF MATTER 

Whilst in Einstein- iEther we can exhibit the large class 
of solutions of the vector field in (fTT)) which lead to caus- 
tic singularities, it is far from obvious that the same oc- 




FIG. 1: Flat Space evolution of a Gaussian perturbation to the 
radial component of the vector field for the Einstein-Maxwell 
case (c2 = 0, no divergence term included). The solution 
displays caustic instabilities. 
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curs in TeVeS. For these solutions, F^. 



and this leads 



geometrically to the integral curves of the vector simply 



FIG. 2: Flat Space Same as figure [T] but for C2 = 0.2 i.e. a 
divergence term included in the kinetic part of the Einstein- 
jEther vector action. The addition of a divergence term sup- 
presses the formation of caustics. 



being geodesies that fall into the gravitational wells cre- 
ated by the matter. However in TeVeS there is direct 
coupling between the vector and matter, and hence in the 
presence of matter one can not have F^^, = 0. If F^^^u ^ 0, 
then the curves of the vector do not follow geodesies, and 
we cannot argue that they must cross forming caustic sin- 
gularities, even if we suspect they might. 

However we can make some analytic progress in the 
absence of matter so T^^, = 0. Then we may consis- 
tently truncate to solutions with constant scalar field. 
Initial data with constant scalar and vanishing scalar 
time derivative, i.e. = on an initial Cauchy surface, 
evolves to have constant scalar. 

We stated earlier that we are interested in the strong 
acceleration regime, for example we look at the dynamics 
on solar system scales or smaller. We denote this scale of 
interest by L. Since the acceleration regime is actually 
determined by the scalar gradient, a precisely constant 
scalar corresponds to exactly the opposite, the low accel- 
eration MOND regime, even if all other dynamical fields 
have characteristic scales given by L. However, we are 
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envisaging a physical situation in which the TeVeS scalar 
has long wavelength fluctuations set by surrounding mat- 
ter, for example set by the galaxy in which the region of 
interest is embedded. These fluctuations will be taken to 
have gradient large enough to place our region of interest 
into the Newtonian regime, which translates to the con- 
dition that the scalar should vary on lengths set by the 
TeVeS scale In Appendix [XI we show that given this 
vast separation of scales, L <c ^, a solution to the TeVeS 
equations where we instead take exactly constant scalar, 
ignore T and set /j, = 1, is a good approximation to the 
full TeVeS equations, within the scale of interest L. [g^l 
For a constant scalar the vector equation reduces to 



gravity wave spacetimes and black hole exteriors (with 
constant scalar). 

We now briefly review some basic facts in GR. For a 
congruence of timclike geodesies, parametrized by proper 
time T with tangent vector fleld with ^''^^^ = —1, we 
may define a tensor field 



(21) 



which then satisfies B^j^^i^^ ~ — 0. We define the 

expansion 0, shear cr^iy, and twist w^jy as. 



KWpF''^" + XA" = 0. 



(20) 



Consider starting with initial data on a Cauchy surface 
H a.t t — where F^i, — 0. Since A'' is timelike, the t 
component A* cannot vanish, and hence the t component 
of the above vector equation sets A = on this initial 
data surface. However, the remaining components i = 
1 . . . 3 then determine dtF*^ = at t = 0. Furthermore 
the Bianchi identity for F^,y, ^[tiFi^a] — implies that 
dtFij = at t = 0. Together these imply that dtF^„ = 



(22) 



where = g^^, + S^^S^^, is the projector onto the tan- 
gent space orthogonal to the timelike geodesies. Then 
Raychaudhuri's equation is. 



de 

rf7 



a^.a'^" -u^^Lo'^" ~R^,eC- (23) 



and so F,,^ remains zero when evolved off the surface E . Now we consider applying this result to our situation. 
Hence we see that starting with initial data d^cj) = F^^ = Recall that since our solutions have F^^ = 0, then the 

integral curves of are timelike geodesies, and moreover 
the tangent vector A^ has unit norm. Hence we may take 
the above to be A^. Then since F^j^i, = the twist 
bjfj_i, vanishes, and hence the congruence is hypersurface 
orthogonal. Furthermore we have = 0, and using the 
fact that a^ya^'^ > therefore arrive at the expression, 



on E implies (in the absence of matter) that the scalar is 
constant, i^^;^ = and A = for all t. As discussed above, 
for F/^ii, = the vector can then be written as = 
with the timelike constraint A^^A^ = — 1 giving the p.d.e. 
in equation (jlSp which can be used to evolve x- Hence 
the initial data for the vector can be parametrized by 
the function x on S which determines the direction of 
the vector on S. 

The dynamics of TeVeS in this truncation are the same 
as those of Einstein-iEther with Maxwell kinetic term 
and no matter. Hence as we claimed above, for F^i, = 
where vector integral curves are timelike geodesies, we 
should expect to be able to form caustics. While this 
is true, and indeed figure [1] gives an example in 1-1-1 for 
the Minkowski geometry, there is no matter to focus the 
geodesies and hence it isn't obvious how generic caustic 
formation would be. If we start with initial data in this 
class of solution - i.e. suitable data for the metric, to- 
gether with the function x which specifies initial data for 
the vector - are the initial data that develop to a singu- 
larity generic, or a special case? We now address this 
by precisely characterizing when initial data will form a 
caustic. We note that while our analysis is given in the 
context of TeVeS, precisely the same argument can be 
made in the context of Einstein- iEther theory, although 
we know of no previous literature doing so. 

Hence we consider TeVeS in the absence of matter, 
with constant scalar, and with Ff^^ — 0. We note 
that the TeVeS vacuum, with Minkowski geometry and 
A^^ = {dt)^ is in this class, and hence we may regard 
the class as a restricted (although not necessarily small) 
deformation of the TeVeS vacuum. The equation for the 
metric immediately gives i?^^ — 0. Thus the class covers 



de 

d^ 



< 



1 



(24) 



where d = (V^A^)(.g^'' + A^'A•') = V • ^ since A^'A|, = 



— 1. Thus we have the result, 



(V.A)-i(r)>(V.A)oi + -r, 



(25) 



along a geodesic with W ■ A — (V • ^)o at the point where 
the geodesic intersects the initial Cauchy surface E. So 
we conclude that if V • A < anywhere on the initial 
hypersurface a, within a proper time — 3(V • A)q^ , V • A 
diverges, signaling that the geodesic congruence ends at 
a caustic singularity. 

In summary, we have obtained the following result: In 
the absence of matter smooth initial data with 8^4) — 
Ffii, = on a spacelike hypersurface E will evolve to 
form a caustic singularity ifV ■ A < anywhere on E. 
Note that while these are solutions with exactly constant 
scalar, the timescale of caustic formation is set by L and 
so they will still be good approximations to the TeVeS 
equations in the strong acceleration regime, as discussed 
in Appendix [Xl 

Whilst this class of solutions is clearly restricted it is 
still physically reasonable and in particular includes ini- 
tial data close to the TeVeS vacuum. The initial data 
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includes the initial data for the metric and for the vec- 
tor, the function x on E. The condition that V ■ A < 
at any point on the initial data surface is very weak, 
and certainly generic within our restricted class. For ex- 
ample, consider the small perturbation from the TeVeS 
vacuum where the metric is taken to be Minkowski and 
the vector near the initial surface E at i = is given by 
X = —t + Sx, for small 5x- Then the singularity condition 
on E, (V • A) ~ V?5x < will be generically satisfied in 
the region surrounding a maximum in % on E. 

Hence our result very clearly highlights the fact that 
caustic singularities do indeed occur in TeVeS. Note that 
we have bounded the time to form the singularity by the 
initial data - no dimensional constants enter that might 
be used to tune away the singularity. It is worth empha- 
sizing that if we were considering a gauge theory where 
Af^ were a gauge potential, then such caustic singular- 
ities would not be a concern. One could simply gauge 
them away. However, the Lagrange multiplier and cou- 
pling of the vector to the matter here mean there is no 
such gauge symmetry, and a singularity of the type here 
is a true physical singularity. The classical evolution of 
the theory is ill-posed once the singularity has formed. 

In a sense we are done. Whilst in the absence of mat- 
ter for more general initial data we have no argument to 
show caustics form, it certainly seems likely. This fact 
is already sufficient to render the dynamics of the TeVeS 
theory rather dubious. However since phenomenologi- 
cally it is the interaction with matter that is of tanta- 
mount interest in TeVeS it is still interesting to under- 
stand whether matter ameliorates the situation or instead 
makes it worse and we spend the remainder of the paper 
investigating this. 

An interesting and possibly related topic is the lin- 
ear instability noted by Seifert (3^ when considering 
perturbations of spherically symmetric backgrounds in 
Einstein- iEther and TeVeS using the methods of [ssl ]. 
Clearly the phenomenon of caustic singularity formation 
is essentially non-linear and one would not expect to see it 
in linear theory. However the dynamics of such singular- 
ity formation may be associated with a growing unstable 
linear mode for the vector field, and it would be interest- 
ing to investigate whether there is indeed a link between 
our discussion and Seifert 's linear instability. Seifert has 
discussed modifying TeVeS in the same sense we do later, 
and his work led to Skordis considering the same mod- 
ification we present later in the context of cosmological 
perturbation theory. 



III. VECTOR FIELD DYNAMICS AND NEW 
STATIC BLACK HOLE GEOMETRIES 

In this section we find new classes of black holes in 
TeVeS where the Einstein metric is static, but the vector 
field has non-trivial dynamics, and its integral curves fall 
through the horizon. The previous static black hole ge- 
ometries of Giannios [l3| and of Sagi & Bekenstein [23 | 



have had the vector aligned with the killing time, which 
naively suggests the vector field in the exterior regions to 
slowly moving matter might dynamically wish to align 
with the matter's natural frame. However, the existence 
of our new solutions clearly shows that this is not to be 
expected, and that the dynamics of the vector exterior 
to a region containing matter may be very complicated, 
and in particular is likely to want to fall towards the 
matter and may form caustic singularities. For an actual 
matter source rather than a black hole where the integral 
curves cannot disappear through a horizon, accumulation 
of curves as the vector falls towards the matter are more 
likely to form singularities. Indeed in later sections our 
simulations with matter show that this is the case. 

Later in the paper we will suggest a modification of 
the TeVeS theory to avoid caustic singularity formation. 
It is worth noting that the black hole solutions presented 
below will not be solutions in this modified theory and 
hence we avoid going into detailed phenomenology for 
these solution here. It would be interesting to study black 
hole solutions in the modified theory we suggest, and we 
make some comments on this in the concluding discussion 
section I VII 

As discussed in the previous section we will be inter- 
ested in length scales L <^ £, where the TeVeS scalar 
is approximately constant in the region of interest but 
varies enough to place the region in the strong accelera- 
tion regime. Thus, as discussed in Appendix [XI we may 
consider solutions to the TeVeS equations with exactly 
constant scalar, ignoring and setting /it = 1 as a good 
approximation. 



A. New static black hole geometries with constant 
scalar 

It was shown by Giannios in fTT*! that if the vector field 
A is aligned with the time translation Killing vector, then 
the solution for is singular unless it is constant. Letting 
that constant be (j)c, the solutions take the form. 



A" 



1 



,0,0,0 



(26) 
(27) 



ds^ = -T{r)df + R{r)dr^ + {dO^ + sin^ Odip^) , 

(28) 

For our new solutions we again have a static Einstein 
metric and constant scalar, but now take A^ ^ 0, 



A'^ ^{A\r),A^r), 0,0) 



(29) 



The scalar field equation is trivially satisfied for a vacuum 
spacetime. The vector field equation becomes 

+ AA" = 0.. (30) 

For the case where a — r the first term vanishes leaving 

AA'' = 0, (31) 
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and thus for A'" ^ we have A = 0. In this case the field 
equations become those of Einstein-Maxwell theory for a 
particular choice of gauge. Given this we expect to find 



Reissner-Nordstrom (RN) black holes, and one can check 
this is indeed t he g eneral solution - we give the argument 
in appendix |B] [63| . We find, 



Diag 



A,, = 



2M Q"" 



<^ = -t± dr\l 



2M Q 



2M Q2 



r^r^ sin^ (61) 




(32) 
(33) 
(34) 




dinates 



FIG. 3: Penrose Diagram for the Schwarzschild MF Spacetime 
where the vector field is aligned with the Killing time as in 
the Giannios solutions IITII. 



and hence see that the solution is indeed simply RN 
where the gauge freedom for the field A, specified by 
the function $, has been fixed up to a sign by the La- 
grange multiplier constraint that A has unit timelike 
norm. Note that for Q — this gauge transforma- 
tion $ is the Lemaitre time coordinate, i.e. the coor- 
dinate time experienced by in- falling geodesic observers. 
Clearly, Ff^^, = when Q = and so the integral curves 
of A are geodesies of the EF geometry. 

One may show that the matter frame metric is RN too, 
after an appropriate coordinate transformation. We cast 
g into Schwarzschild-type coordinates, (r, p) defined as, 

t{r,p)^e-^^T + f{p) r{p)=e^^p, (35) 



where / is given by. 



gtt ' 



(36) 



These transformations leave all components r indepen- 
dent, and g in standard RN form in Schwarzschild coor- 




where now the mass and charge are given as. 



= M = e-^'f'-M - (1 




j^Q, (38) 



+ (1 - e-^^=) ) Q2. (39) 



This shows the metric is still RN, but with shifted horizon 
positions. For example, in the Schwarzschild case {Q = 
0) we have that the MF and EF horizons are related by 



(40) 



Bekenstein |8| demonstrates the speed of scalar pertur- 
bations (at fixed background vector) and vector pertur- 
bations (at fixed background scalar) is sub-luminal with 
respect to electromagnetic propagation only if > ev- 
erywhere. So we conclude that in this case, the horizon 
seen by the matter fields is smaller than the horizon for 
the gravitational fields (0, A, g). 

We now approach TeVeS solutions from another direc- 
tion. Starting with the standard RN solution in Einstein- 
Maxwell theory, we may obtain a vacuum solution of 
TeVeS with constant scalar and A = provided we may 
choose a gauge such that the vector potential satisfies 
the TeVeS constraint A'^A^ = — 1. Hence we find a large 
class of solutions with metric ([5^ and vector of the form 
(1551) but where, as in the earlier equation (fT5|) . we regard 
the gauge condition as a first order p.d.e. in time for 



2 Q 



1 



K r ' (-.9**) 
V(-.9")(l + 9'*9^$) + (.g"a,$)2 



(41) 



Hence we may take the solutions to be characterized by 
the charges AI, Q and also ^{t = 0), from which we evolve 
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FIG. 4: The Penrose Diagram for the generahzed 
Schwarzschild MF Spacetime (Q = case) with the EF hori- 
zon position. In this case the vector field is free-falling along 
geodesies of g through both horizons. 



to construct $ for all t. The solution above in equation 
([M)) is a stationary solution to this p.d.e. for $. In 
general however the solutions to this p.d.e. will have 
complicated time dependence. 

Note that for Q — 0, this class of solutions has = 
F^u = as well as constant scalar and hence falls into 
our earlier class discussed in section|lIl Thus in this black 
hole background we can again precisely characterize caus- 
tic singularity formation by the previous statement that, 
letting (V • A)q be the value of V • A on a surface E, then 
if (V • A)o < at any point on S a singularity will form 
within a proper time — 3(V • A)^^ measured along the fu- 
ture of the integral curve of A through that point. Since 
this condition is generic within our class of solutions for 
Q = 0, we see that caustic singularity formation is to be 
expected in the exterior of these black holes. Whilst we 
have no argument that the same is true for Q 7^ 0, we 
expect it is likely. 

The presence of the matter appears to attract the vec- 
tor field integral curves. Naively this focusing would 
seemingly make caustic singularity formation more likely. 
However, interestingly the presence of the horizon actu- 
ally renders the singularity formation less severe in the 
sense that if the time to singularity formation is sufficient 
that the integral curve of A has fallen inside the horizon, 
an external observer need not care. Indeed in the sta- 
tionary solution (|34p above this precisely happens, with 
a caustic singularity occurring at the black hole singular- 
ity itself. Of course, if the matter source is not a black 
hole, but rather a compact object without horizon, then 
we still expect the vector field curve attraction, but now 
there is nowhere for the curves to go, and hence the ex- 
pectation that the matter focuses the vector to form sin- 
gularities would hold. Later in the paper we investigate 
this. 

We conclude this section by commenting that we have 
examined the case of constant scalar, stationary black 
hole solutions. One might wonder whether stationary 



solutions with non-constant scalar can be found, which 
share the symmetry and asymptotics of those found here. 
We address this questions in the Appendix [Cl finding ev- 
idence that no solutions exist near to the ones above with 
constant scalar. We show that for a linear perturbation 
of the scalar about our constant scalar solution a singular 
develops exterior to the horizon, and that performing a 
full non-linear numerical evaluation of such a solution one 
finds both the EF and MF metrics are nakedly singular. 
We cannot argue that no non-constant scalar black holes 
exist [g^l but do expect there are none that are quali- 
tatively similar to the constant scalar solutions we have 
found. Note that this is compatable with the argument 
presented in Appendix [X] since the kind of non-trivial 
scalar considered there need not be static nor share the 
same asymptotics or symmetry. 



IV. VECTOR FIELD DYNAMICS AND 
MATTER: NUMERICAL SIMULATION 

In section [ll] we have shown that in the absence of 
matter, a large class of deformations of the TeVeS vac- 
uum initial data quickly terminate in caustic singularities 
upon time evolution. In section IIIII we have shown that 
black hole solutions of TeVeS have complicated vector 
dynamics, which again include caustic singularities, and 
in particular that the black hole appears to attract the 
integral curves of the vector field toward it. We might 
then expect this to occur for any matter source, and then 
imagine that such an attraction which focuses the inte- 
gral curves is likely to generate caustic singularities. This 
is too quick however, as matter couples to the vector field 
in TeVeS in a complicated fashion, and hence we have lit- 
tle intuition or analytic control over what happens. 

It is the subject of this section to investigate the vec- 
tor dynamics in the presence of matter using full nu- 
merical evolution of the TeVeS equations of motion. To 
make progress we restrict ourselves to spherical symme- 
try. We are then able to consider both gravitational col- 
lapse of a matter scalar field, and evolution of an initially 
near static boson star. In both cases we find the vector 
curves in the region exterior to the matter are indeed 
attracted towards the matter and do form caustic singu- 
larities. One might be concerned that imposing spherical 
symmetry restricts to a rather special class of solutions 
which focuses energy at the origin of spherical symmetry. 
However the caustic singularities actually form away from 
the origin and hence the singularities themselves locally 
have a planar symmetry, and seem not to result from the 
peculiarities of spherical symmetry. In both cases the 
TeVeS scalar is fully dynamical and non-constant, hence 
justifying the Newtonian regime approximation (/i — 1, 
neglecting J-). It remains smooth where the caustics form 
indicating that the scalar plays no role in the pathalogical 
vector dynamics. 

Full details of the numerical implementation and con- 
vergence and constraint tests are postponed to Appen- 
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dices [D] and m 

A. Scalar field collapse 

Our first matter system is the collapse of a complex 
scalar field. We perform an integration of the field equa- 
tions from an initial spherical shell of scalar matter. We 
use a canonical complex scalar field Xi whose action is 
constructed using the matter frame metric g. Using time 
symmetric initial data and a radial gaussian shell for the 
matter field x, the matter energy density will split into 
ingoing and outgoing components and for sufficient am- 
plitude of the initial shell, the ingoing component will be 
focused at the origin into a high enough energy density 
to form a black hole. We use the coordinate system 

ds'' = -T\t, r)dt^ + e«(*-'-)dr2 + r^dVil (42) 

which clearly covers only the exterior region of any black 
hole that might form. We emphasize that since the mat- 
ter X couples to the TeVeS vector and scalar, both must 
be evolved and have non-trivial dynamics. 

For the metric we choose time symmetric initial data 
that satisfies the constraint equations. The TeVeS scalar 
we take to be constant initially, and time symmetric, and 
for the vector, we take — a.t t = and choose 
such that F^i^ = at t = 0. However, since there 
is matter, in contrast to the discussion in the previous 
sections above, F^^ will immediately evolve to be non- 
zero, and </) to be non-constant. 

We find that for different initial Gaussian shells of x, of 
sufhcient amplitude to ensure non-linear dynamics when 
the ingoing pulse reaches the origin (for weak amplitudes 
the energy density simply passes through the origin and 
radiates to infinity as is the case in usual standard gravity 
coupled to a scalar) rather similar qualitative behaviour 
results. Figures [5l Eland [7] give the results of a representa- 
tive evolution. Even though i^^jy 0, we do see that the 
vector field curves are initially attracted to the matter 
shell. Evolution proceeds with the ingoing pulse deform- 
ing the geometry as for a usual scalar collapse. However 
before a horizon can form - recall our coordinate system 
only covers the black hole exterior - we see the formation 
of a caustic singularity. This is signaled by the divergence 
of A which develops a growing spike on constant t slices 
as seen in figure [S] At this point, dynamical evolution is 
no longer well defined. As noted above, the singularity 
forms well away from the origin of spherical symmetry. 
Figure [7] shows that the TeVeS scalar remains small and 
smooth up to this point, which suggests that dynamically 
it does very little, if anything, to prevent singularity for- 
mation. It is interesting that whilst we have found many 
candidate black hole end states for such a collapse in the 
previous section, the actual dynamics of the collapse is 
sufficiently badly behaved that we cannot even see an 
apparent horizon form. 

One might worry that the choice of initial data with 
Ffii, = is somewhat special (even though evolves 



to be non-zero). [65| Indeed, since our earlier analytic 
arguments were for F^i, = it is useful to check that 
caustics may also form for initial data with ^ 0. 
Another question is whether the magnitude of if, n play 
a role in the singularity formation. For these reasons we 
present the result of another simulation in figures [SJ [5] 
and [ini These simulations used initial data for the vec- 
tor where Ar = Ar — 0. Hence the initial data as a 
whole is now time symmetric, and ^ initially. For 
the simulation shown we have also taken larger K, k. We 
observe caustic formation again. Indeed the singularity 
forms earlier. Experimentally we find that for larger K, n 
a caustic forms earlier, which is to be expected as the vec- 
tor is more strongly coupled to the dynamics of the other 
fields. Thus we see that for two very different choices of 
the vector initial data caustic singularities result. 

B. Perturbations to a boson star 

Scalar field collapse is an extreme dynamical process 
which is highly relativistic. It is interesting to consider 
whether an initially non-relativistic matter source also 
seeds an attraction of vector curves and subsequent caus- 
tic singularity. To this end we examine the full dynamics 
of TeVeS in the presence of an initially quasi-static boson 
star [34] . 

We use an identical numerical method and boundary 
conditions to integrate the field equations as for the scalar 
collapse above. To create the boson star we follow Gleiser 
(35| . as a matter source we use a complex scalar, x, now 
with potential T^dxD- We begin by finding a static boson 
star solution. This is achieved by imposing the following 
separable solution to the x equation of motion with the 
potential l^(|xl) = '^t^xx 

x{t,r)=xo{r)e''^\ (43) 

where xo is real. While we term this a 'static' star, we 
note that in fact x has the above phase rotation, although 
all other fields are indeed static. The metric functions 
T, R are taken to depend only on r, and likewise for the 
TeVeS scalar. The TeVeS vector is chosen to be aligned 
with Killing time, so = (-r(r), 0, 0, 0). The radial 
profile for each of these functions is obtained via a shoot- 
ing method - we fix the value of Xo(0) a-nd m so that the 
resulting solution will have flattened out well before the 
boundary of our numerical grid. We then fine tune the 
value of w to obtain the profile for xo with no nodes; this 
is the ground state star. We also choose the parameters 
so that the start has a low density compared to its radius 
and hence T{0,r) ~ 1, so the backreaction of the star is 
weak - it is non-relativistic. Note that for this static star 
Ffiu 7^ 0, and the TeVeS scalar is non-constant. 

To consider a dynamical perturbation of this static star 
we take similar initial data to the scalar collapse. We 
take T, R, the TeVeS scalar and the matter scalar to have 
initial data simply given by that on a constant time slice 
of the static boson star solution above. However we now 
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FIG. 5: Massless scalar field shell-collapse in an attempt to form a black hole. Upper Frame: The light cone structure is shown 
overlaying contours for the amplitude of the x matter field. Lower Frame: Integral curves for the TeVeS vector A. In this 
evolution we see a caustic singularity develop around r ~ 11 at f ~ 3.2. The light cone structure indicates an apparent horizon 
is likely to have formed if the evolution had not been terminated by the caustic. The initial data for the vector was Ar = 
with Ar chosen so -F^i/ = 0. The coordinate system is (|42)). with parameters K = k = 0.01 



take A"" = and F^^ = at t = (although again 
F^i, will not evolve to remain zero due to the boson star 
matter) . Thus the vector is not now aligned with Killing 
time and dynamics will ensue. Note however that since 
the stars considered are in the non-relativistic regime, 
this perturbation to the vector field initial data is small. 
Hence the evolution is a non-relativistic process in its 
early stages . 

We performed evolutions for a variety of star config- 
urations, obtaining qualitatively similar results. A rep- 
resentative choice is illustrated in figure [TTl [12] and [131 
This shows that for our quasi-static initial configuration, 
the vector field curves fall in towards the star and do 
evolve to form a caustic singularity. However the sin- 
gularity does not form in the interior of the star as one 
might naively expect. This is essentially due to the vec- 
tor coupling to the matter, which apparently leads to a 
repulsive effect as we see the integral curves are clearly 
repelled from the origin of the spherical spatial geometry. 
While F^^ is not zero outside the star, as it is sourced 
by the stellar matter, and then the region where F^i, is 
non-zero propagates outward, we see that it does not stop 



the vector curves falling towards the star and eventually 
'colliding' with the curves that were 'bounced' out of the 
interior of the star. The singularity appears rather sim- 
ilar in nature to that in the case of the scalar collapse, 
and figure [12] clearly shows a growing vector divergence 
at a finite radius as we approach the caustic. Figure [T51 
shows that once again the TeVeS scalar remains small 
and smooth up to the singularity. 

Thus we have seen in this section that even start- 
ing with initial data whose short term evolution is non- 
relativistic, pathological vector behaviour may quickly 
follow. In particular we see visually that since the in- 
tegral curves are, at least initially, following an approxi- 
mate timelike geodesic, the timescale for this singularity 
formation is of order the gravitational in-fall time. Thus 
in a Newtonian, quasi-static regime such as the neighbor- 
hood of the Earth, one might expect caustic singularity 
formation to occur on the order of hours, and after that 
point classical evolution is ill defined. This poor dynam- 
ical behaviour is clearly a serious obstruction to consid- 
ering TeVeS as a phenomenological theory of modified 
gravity. 
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FIG. 6: —V • A for the simulation of massless scalar field shell-collapse, as in figure [5] Spiking of the 4-divergence of —A 
corresponding to the convergence of its integral curves is easily seen, signalling the formation of a caustic. 



V. MODIFYING TEVES TO GET A WELL 
BEHAVED VECTOR DYNAMICS AND A MOND 
LIMIT 

Having demonstrated the formation of caustic singu- 
larities in TeVeS in various contexts which render the 
classical dynamics of the theory quickly ill defined, we 
now propose a correction to the vector part of the action 
dH) which may ameliorate this problem. The problem 



J 



is essentially due to the Maxwell structure of the vec- 
tor action. There is no energy cost when the divergence 
of the vector becomes large. Our modification is sim- 
ply to introduce terms that disfavour large divergences. 
We simply take the vector action to be the most gen- 
eral diffeomorphism invariant action which is quadratic in 
derivatives and consistent with the = —1 constraint. 
This action is of course precisely the one used for the 
vector in JE theory. We begin with 



167rG 



C2(V^A'^)2 - c^A^At^ + XiA^ + 1) 
I 



{-gy/^d\ (44) 



where c+ = ci + cs, a - cs = 2K, A^^ = A^V^Af^ and 
Sfj.u = V^Ajy -I- Vi/Ap. We retain the pure TeVeS scalar 
action. The metric redefinitions which may performed 
in JE theory [361 to remove one of these terms are no 
longer applicable here, as there are no non-trivial field 
redefinitions which leave the scalar action form-invariant 
at the same time as retaining the unit-norm constraint. 
As noted earlier, Seifert [s^l has already proposed such a 
modification of TcVcS motivated by finding a linear insta- 



bility about certain spherically symmetric backgrounds 
(those of Giannios QJ^) and leading from this Skordis [1^ 
has recently derived the equations of motion and studied 
the cosmological perturbation equations. Our emphasis 
here is to check that the MoND limit is still recovered 
with this modification, and this has not previously been 
addressed - without this, of course, the modified TeVeS 
would be unlikely to provide any alternative explanations 
for dark matter. 



13 




-0.Q0990 



-O.OIOQO 



FIG. 7: —(j) for the simulation of massless scalar field shell-collapse, as in figures [5] and |6l The scalar field remains small and 
smooth up to the formation of the caustic. 



With this modified vector action, only the Q term of 
the metric field equation is affected, and we obtain [66l | 



J 



+C2 (g^..V^ (A^V • A) ~ A(,,V,)V -A^^iV-A) 



A A^Ay 



-A„A' 



f C4 (i^i^ + A^A^f,Vy)A'' - V„ 
KVf.F^"' + y V^S*^"^ + caV (V • A) - CiA-^VA^ + aV^ (i'-A'^) 

I 



(45) 



for the vector and scalar equations. We will refer to 
Bekenstein's theory as 'pure TeVeS', and the theory with 
the modification as 'modified TeVeS'. 

The new parameters introduced into the action, 
c+, C2, C4 will certainly be subject to physical constraints. 
The considerations are likely to be similar to those con- 
straining the Einstein- jEther parameters, reviewed re- 



cently by Jacobson 31]. One complication is that the 
physical fluctuation modes of the vector now have dif- 
ferent wave speeds when the more general vector ac- 
tion is introduced. In particular this leads to new con- 
straints from Cherenkov radiation produced by cosmic 
rays [37, 38], although since the matter couplings are dif- 
ferent from those of Einstein- jEther such analysis would 
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FIG. 8: As for figure O although now with initial data Ar = Ar = 0, so that F^i, initially. This simulation was performed 
with larger K = k = 0.1. We see again caustic formation, now much sooner. 



likely have to be repeated for the modified TeVeS. We 
leave determination of constraints on these parameters 
for future work. 

We have modified the TeVeS vector action in the hope 
that it will alleviate the problem of singularity formation. 
An absence of caustic formation is a crucial requirement 
for the theory to be dynamically well behaved, though we 
do not attempt to assess whether this is actually the case 
for our modified theory. We leave this as an interesting 
open problem. Now, assuming that the dynamics of this 
theory are in fact good, we would then require the the- 
ory to have the appropriate phenomenology. That is, we 
would like modified TeVeS to have inherited pure TeVeS's 
MoND limit. In order to check whether this is the case, 
we perform a Newtonian analysis where the goal is to 
obtain the equivalent of Poisson's equation for MF New- 
tonian potential. We perturb the EF metric to leading 
order in the Newtonian expansion as. 



9tt 



-1 - 2V 



5y = (1 - 2V)6^j gu = 0, (46) 



where we linearize the equations in V, and ignore time 
derivatives at leading order. The matter source has only 



the non-trivial component Ttt — p. As for the EF metric, 
for the MF metric we take, 



~gu = -1 - 2$ = {l- 2$)5, 



gu 



0, (47) 



and in the Newtonian expansion only the time compo- 
nent of A is non-trivial at leading order, and is deter- 
mined from V by the condition A^'A^ = -1. The TeVeS 
scalar is written as (j) = (pc + S(f>. In the Newtonian limit 
at leading order we take Scj) << 1, and again neglect 
time derivatives. One may then verify that the disformal 
relation then relates these perturbations: ^ — V + cj). 

Consider first the scalar field equation at leading order 
in the Newtonian expansion. 



(V<5<^)') 



nGp. 



(48) 



Note that we have not linearised the argument of the fi 
function in 5(j) since while d(j) « 1, the TeVeS parameter 
I is precisely large enough to balance this - hence we may 
recover MoND at leading order in the Newtonian analy- 
sis. There is no such subtlety in the other field equations, 
and we may straightforwardly perform a linearization in 
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FIG. 9: "V ■ A for the simulation in figure |51 



the Newtonian potentials. All components of the vector 
field equation vanish except for the f-component, 



A + - ^) W = (-4</.,)87rGp, 



(49) 



and similarly all components of Einstein's equation van- 
ish with the exception of the tt-component, 



A+ (2 + C4 -c+)VV = (1 



inGp. 



(50) 



where we note that, as shown by Bekenstein [8|], the term 
doesn't contribute to the stress tensor in the leading 
order Newtonian limit. Combining the above results, we 
learn how $ is related to <I>Ar and 4>, 



i + ^^^^:^-40,)$A, + 0. (51) 



Hence, by following the same arguments that apply to 
pure TeVeS theory (as discussed by Bekenstein in [s*]), 
MoND phenomenology results from and ([5T|) when 
C4 and c+ are suitably constrained. Hence our proposed 
modification of TeVeS, which likely can avoid caustic for- 
mation for specific parameter ranges, does indeed cor- 
rectly reproduce the MoND limit which is the raison 
d'etre for TeVeS. 



VI. SUMMARY AND DISCUSSION 

We have argued that Bekenstein's original formulation 
of TeVeS, while reproducing MoND phenomenology, is 



actually dynamically badly behaved. We have shown an- 
alytically and numerically in a variety of situations that 
the integral curves of the vector field quickly and gener- 
ically evolve from regular initial data to caustic singu- 
larities. Once this occurs the classical evolution to the 
future is ill-posed. 

Since the time scale to a singularity is not suppressed 
by any parameters, and is only determined by the initial 
data itself, it seems the original formulation of TeVeS is 
unlikely to provide a realistic theory of modified gravity. 
Put another way, while TeVeS does reproduce Newtonian 
and Modified Newtonian dynamics in a non-relativistic 
regime, it appears that in many cases this regime is un- 
stable, the instability leading to the caustic singularities. 
We stress that already the analytic arguments of section 
im whilst made in the absence of matter, already high- 
light instability in the dynamics of TeVeS. The latter 
sections of the paper merely serve to illustrate that the 
situation remains unchanged when one considers the dy- 
namics in the presence of matter. Indeed since matter 
can focus vector field curves towards it, it can make the 
situation worse. 

It is useful to contrast this situation with the singular- 
ities that form in GR. We are very familiar with the fact 
that given certain initial data, matter can collapse to a 
singularity in GR, on the time scale of the gravitational 
in-fall time. However, in GR Cosmic Censorship means 
that these singularities are always hidden behind event 
horizons. Hence evolution outside the horizon is perfectly 
well defined. In contrast, the vector field singularities we 
have exhibited here lie outside any event horizon, and 
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FIG. 10; —4> for the simulation in figures [8] and |9l 



hence evolution is impossible in the future null cone of 
these points unless there is some way to understand these 
singularities beyond the TeVeS effectiye field theory. 

A similar situation arises in the perfect fluid, dust 
description of dark matter. The eyolution of such a 
fluid also forms caustics (shocks) on time-scales of the 
in-fall time. Caustic formation in this case signals the 
breakdown of the fluid description of the dark matter 
and a requirement for a microscopic, particle description. 
Similarly our results indicate that the dynamical regime 
where TeVeS can be applied is limited by the breakdown 
of the effectiye theory on in-fall timescales. In this case 
howeyer we don't haye a microscopic description to tran- 
sition to and the regime where the effectiye description 
appears to break down is releyant to the motiyation of 
the theory itself. 

Nonetheless, we believe a relatively minor modiflca- 
tion, namely generalizing the vector action of TeVeS to 
a form like that of Einstein- jEther is likely to be able to 
give a dynamically well defined theory, which as we have 
shown still gives MoND behaviour for non-relativistic sit- 
uations. It is interesting that a possibly related instabil- 
ity was observed for linear perturbations about the spher- 
ically s ymm etric static backgrounds of Giannios |17| by 
Seifert [32| and the same modification was proposed, al- 
though the recovery of the MoND limit had not previ- 
ously been checked. 



We emphasize that the detailed predictions of this 
modified TeVeS will likely differ from the original TeVeS 
theory, and therefore any phenomenological studies of 
TeVeS testing its ability to explain astrophysical or cos- 
mological data without dark matter should be careful to 
include the necessary modification. It would be interest 
ing to revisit the questions addressed in 
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We have given large classes of new black hole solutions 
in the original TeVeS theory. However, these and the ear- 
lier solutions of Giannios are not solutions of the modi- 
fied TeVeS theory. Instead in the case of modified TeVeS 
black hole solutions with a constant scalar (and it is likely 
there are not 'nearby' solutions with non-constant scalar) 
will be identical to those in Einstein- iEther theory dis- 
cussed by Eling & Jacobson [H^ [6^ . In particular there 
is no 'charge' parameter Q, with the static black hole 
geometries only being parametrized by one parameter, 
the mass M. However the qualitative behaviour of the 
solutions is somewhat similar to ours, with the integral 
curves of the vector field falling into the future horizon. 
It would be very interesting to study the geometry of 
the MF metric for these solutions. As a passing com- 
ment Eling et al. [6l| have argued that such black holes 
may violate the Generalised Second Law. It is interest- 
ing that while the black hole solutions we found in the 
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FIG. 11: A complex scalar field x (grey shading) initially in the lowest static mode of a quadratic potential, as a matter source 
for a full TeVeS evolution. Shown also are the integral curves for the TeVeS vector A, illustrating the phenomenon of caustic 
formation where the initial data is quasi-static. 



original TeVeS do not necessarily violate the law since 
their entropy depends on multiple charges M and Q, for 
black hole solutions in the modified TeVeS there is only 
a one parameter family and the arguments of Eling et al. 
apply. 



Acknowledgements 



be of order planetary or solar system scales. We claimed 
that to any solution with constant scalar, one can con- 
sider deforming the solution by adding scalar gradients 
that are tiny compared to our scale of interest L, and 
hence effectively negligible, but that would still be large 
enough over the region of interest to put the theory into 
the strong acceleration regime. The equation determin- 
ing a is. 



We wish to acknowledge Pedro Ferreira, Andrew Jaffe, 
and Joao Magueijo for useful discussions. We also thank 
Jacob Bekenstein and Martin Feix for valuable com- 
ments. TW is supported by an STFC advanced fel- 
lowship and a Halliday award. BW is supported by an 
STFC studentship. CRC acknowledges support from the 
Nuffield Foundation. 



APPENDIX A: APPROXIMATELY CONSTANT 
SCALAR SOLUTIONS IN THE NEWTONIAN 
REGIME 

In the main text we have made the claim that in the 
absence of matter we may formally consider solutions of 
TeVeS with constant scalar, and yet for the length scale 
of interest, let us call it L, which is much shorter than 
the TeVeS length scale £, we are still in the Newtonian 
regime, /j, ~ 1. We envisage our scales of interest L to 



where we have written /i — kGo^ and y — Ki'^\(j)\^ with 



9 



Q/3 



A'^A'^ and we have introduced the notation 



= h°'^(j>^a4',p- Let us also choose, following Beken- 
stein Q, an ^ such that in the Newtonian limit /x ^ 1, 
we have, 



-ln(l 
2 ^ 



3 1 



2^ ^ 7 



4 1 



and hence in this limit we have, 



(A2) 



(A3) 



We see a potential dilema in this claim is that naively 
/i ~ 1, ?/ — !■ oo appears precisely at odds with a constant 
scalar which has y = 0. Hence how can a constant scalar 
solution ever be 'close' to a solution in the strong acceler- 
ation regime /i ~ 1 . The resolution is that what matters 
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FIG. 12: —V • A for the simulation of a complex scalar field as a ground state Boson star, as in figure 1111 Spiking of the 
4-divergence of —A corresponding to the convergence of its integral curves is easily seen, signalling the formation of a caustic 
singularity. 



is not but and hence one can have a tiny 

gradient on scales of interest L over a region of size L, 
but provided £ is large enough, one can still have /i ~ 1. 
We will now formalize this more carefully, by providing a 
prescription to take a solution of Einstein-Maxwell the- 
ory and then generating a solution of TeVeS with almost 
constant scalar in a controlled manner. 

In the absence of matter the scalar equation reduces 
to. 



0, 



(A4) 



and the vector, which satisfyies the constraint A^A^ ~ 
-1, obeys, 

K^IpF^"^ AA" + —A^(l,^pg°'^(j)^^ = 0, (A5) 

with the Einstein equations governing the metric becom- 
ing, 



Gap = 0, 



scalar 

Q/3 



e 



Q/3 



(A6) 



with, 



e 



Q/3 



EE SttGct^ 



5,Q(p,/3 



2 

2 ^ 



',v 5q/3- 



\u gap 



(A7) 



eif"- = K { F^f'F, 



gapF^ - XAaAp (AS) 



Q^o = -27rG^£-^a^T{KGa^)ga 



/3 



(A9) 



We consider the dimensionless constants K, k to be small 
but order one in what follows. 

Let us now consider a solution, Aa,gap of Einstein- 
Maxwell theory with gauge constraint A^^A^ = — 1 im- 
posed using the Lagrange multiplier A as for TeVeS, 



Gafi=K[ F^^Fp^ - ^ijapF' ) - XAaAp 



K\7gF^'^ + Ai" = 0. 



(AlO) 



Here quantities with hats are composed from Aa^gap- 
Let us characterize the length scales of interest by the 
length L, and restrict our attention to a spacetime region 
V of spatial size ^ L. Hence the curvatures of interest 
in the solution {Aa,gap) will scale as 1/L^ are we are 
not interested in much smaller curvature scales. We are 
envisaging that this scale is of order the solar system or 
less and compared to the TeVeS lengthscale £ we have a 
vast separation of scales, 



L « I. 



(All) 
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FIG. 13: ~<j) for the simulation of a complex scalar field as a ground state Boson star, as in figures [TTI and [12] 



We begin by constructing a solution for a scalar on 
the fixed solution {Aa^gap) in the spacetime region V, 
of spatial scale ~ L. We take the scalar to obey the 
equation, 

yi3{h^''^,c) = (A12) 

where h""^ = g°'^ — A^A^ , and we emphasize that we 
are ignoring any backreaction - it is simply a scalar on 
our fixed solution. We require that in our region V the 
solution everywhere obeys the condition, 

1 » /l"''0,a0,/3 » (A13) 

For example, taking the trivial flat space solution gap — 
rjap, A — -g^, one might choose the scalar to be = ax 
where x is one of the spatial coordinates, and a is a 
constant in the range 1 >> L^a^ >> L/£. In general we 
expect to be able to find solutions obeying the condition 
(1X13)) . although we do not have a formal existance proof 
of this. 

From this Einstein-Maxwell solution and the associ- 
ated scalar solution {(f>, Aa, gap) we may construct an ap- 
proximate solution of the TeVeS equations, {(j>, Aa, gap), 
perturbatively in the dimensionless constant. 



e = 




(A14) 



as, 

A^ = 4 + e4^^ 

A = A + e^AW 

= <t,Q + e 4, + €^ (t>^^\ (A15) 

where is a constant and plays no role in the vacuum 
TeVeS equations which only involve (p derivatives. Tak- 
ing the limit e — s- 0, ie. looking at small scales compared 
to the TeVeS length scale I, we therefore see that the dy- 
namics of TeVeS on scales '--^ L is given by precisely the 
Einstein-Maxwell solution {Aa^gap) ^^nd hence by an ef- 
fectively (although not precisely) constant TeVeS scalar. 
We may think of the e — > limit as fixing L and taking t 
to infinity, or alternatively and more physically for fixed 
^, focussing our interest on smaller and smaller length 
scales L. 

Let us now check this claim. Firstly let us consider 
the scalar equation at leading order in the e expansion. 
Consider the behaviour of From our condition 

(|A13[) above we see that, 

1 » ^ + 0{e^)) » J. (A16) 

so that for e we have, 

f h'"'(l>,a(l),p » \ (A17) 
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and thus we see that, 

fi{Kp\^f) ^ 1 - 0{e^) (A18) 

so that in our region V we are forced into the Newtonian 
regime of TeVeS, even though the scalar field gradient is 
perturbatively small compared to the scale L of interest 
in our region. Hence the TeVeS scalar equation becomes, 

y^ih^^^c^^a) = O(e'), (A19) 

with the right-hand side coming from the non-constant 
part of /i. This is indeed consistent with our ansatz 
(|A15|1 above for the constant and (j> pieces with the cor- 
rection term, (f>^^\ accounting for the lower orders. 0^^^ 
is sourced primarily by the 0(e) corrections to Vq and 
/i"'^, from the metric and vector corrections Aa^ and g^^p, 
with the source from fi actually being sub leading to this. 

We now show that just as the scalar equation is con- 
sistently solved perturbatively in e by our ansatz, the 
Einstein and vector equations are too. In particular we 
must show that the backreaction in the Einstein equa- 
tions from the scalar and TeVeS function are small 
compared to the characteristic curvature scale in 
the solution {g^,^, ^p)- Now following from our condition 
(jAl3P we have that the scalar in our region obeys the 
bound, 

» /l"'3(/.^a</',;3, (A20) 

and in addition we have an estimate for the contribution 
of J- in the Einstein equations, 

~ -^(e^lne) g^p. (A21) 

Hence we see that the Einstein and vector equations to 
lowest order in e reduce simply to the Einstein-Maxwell 
ones, and our ansatz (|A15|) therefore solves them to low- 
est order. The leading higher order corrections come from 
the perturbatively small backreaction of the scalar, and 
lead to the 0(e) corrections to g^^,A^j_, with the TeVeS 
function essentially being negligible in the Newtonian 
limit as discussed in Bekenstein's original paper. 

We have now more carefully justified our claim in the 
main text, namely that we may consider the TeVeS scalar 
to be effectively constant, and still be in the Newtonian 
regime — > 1, provided we are restricting our interest to 
a region of scale L << as we are in the main text. Our 
first application is to consider a bound on caustic forma- 
tion time using Raychaudhuri's equation for the Einstein- 
Maxwell system. Since caustic formation is local, we are 
only concerned with the spacetime in the region of scale 
L where the caustic forms, and not the asymptotic be- 
haviour of our geometry. The physical setting would be 
caustic formation on, for example, solar system scales L, 
with the gradient of the scalar arising from much larger 



galactic scales The second application is to embed the 
Einstein-Maxwell black holes in TeVeS. Again, since we 
are not interested in the far asymptotic region of these 
solutions, we may again employ our approximation to ig- 
nore scalar gradients and the TeVeS function T in the 
stress energy. For a region of size L surrounding the 
black hole, the corrections will be characteristic scale £ 
and for any astrophysical black hole, phenomenologically 
this region of interest certainly obeys L << £. 

We conclude with our previous example; ga/3 — rjap, 
A = A = 0, and scalar solution on this background 
(j) = ax with the constant a obeying 1 >> V^a^ » 
L/£. Hence we may take a = e/L. In this case the 
exact TeVeS solution can be found in the absence of the 
T term, which we have argued is subdominant over the 
other corrections; 

ds^ = -a{xfdt^ + dx^ + -^{dy^ + dz^) 



a{x) dt 

= </'o + e| (A22) 
with the function, 

a{x) = Ae^^V^^^ (A23) 

and one finds A — IGKna^ / {k{2K — 3)). Note that y is 
actually a constant for this solution, y — k e"'''/'* , so that 
fi is also constant. Expanding in e one finds this to be 
consistent with the ansatz. We have ignored the T term 
in this exact solution and a calculation confirms this is of 
order the estimate (jA21[) . and hence vastly subdominant 
to the leading e corrections in the solution. 

APPENDIX B: CONSTANT SCALAR STATIC 
BLACK HOLE DERIVATION 

We begin with a spherically symmetric, stationary sys- 
tem in Schwarzschild coordinates for which, 

= ^c, (Bl) 

= {A\r),A^ir),0,0), (B2) 
ds^ = -T{r)dt^ + R{r)dr'^ + 

[dO^ -f sin^ Odip^) , (B3) 

where 0c is constant. The scalar field equation is triv- 
ially satisfied for a vacuum spacetime. The vector field 
equation becomes 

KVfjF'^"' + AA" = 0.. (B4) 

For the case where a = r in equation (|B4p , the first term 
vanishes, leaving the branch choice 

XA"- = 0, (B5) 
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thus for we have A = everywhere. In this case 

the field equations become those of Einstein-Maxwell the- 
ory for a particular choice of gauge. Given this we expect 
to find Reissner-Nordstrom black holes. Expressing the 



field equations in terms of our metric and vector ansatze, 
we obtain three useful components of the Einstein equa- 
tions, 



(tt) T_^{rR' + R^-R) = K^A',f, (B6) 

{99) (2T'i?r - 2R'T^ - T'R'Tr + T" RTr - T'^Rr) = Kr{A'^f , (B8) 



where primes indicate derivatives with respect to r. Elim- 
inating between equations (jB6|) and (|B7p leads to the 
relation 



TR = Ci, 



(B9) 



where Ci is a constant, which we set to 1 using the free- 
dom available in rescaling the t coordinate at this stage. 
Performing the same elimination between equations (|B7[) 
and (jB8|) and substituting for T using (jB9p , we arrive at 
a solution for the second metric component. 



i?^(l-^ + ? 



with the corresponding solution for T, 



2M Q2 
T= ( 1 + ^ 



(BIO) 



(Bll) 



as r — !■ oo, must be driven to zero so that isotropy 
is restored at spatial infinity. This expression therefore 
forces Af ^ 1, and so we find C2 — ±1. We choose the 
vector A to be future pointing at spatial infinity, and so 
we pick C2 = —1. The value of C3 can then be deter- 
mined straightforwardly by substituting the expressions 
back into the Einstein equations. In particular, for the 
(9, 9) component we find 



KCl 
2r2 ' 



(B15) 



justifying our choice of the positive quantity . We then 
identify 



C3 



(B16) 



where M and are integration constants. The quantity 
will indeed turn out to be positive. 
Consider the t component of the co-vector field equa- 
tion 



1 



{2KrAlRT - KrA^T'R - KrA[R'T+ 

iKA^RT) = 0. (B12) 



Note that the radial vector component appears nowhere 
in these field equations, and will be determined alge- 
braically using the field equation for the Lagrange mul- 
tiplier field, A. Substituting in the metric components, 
and solving the resulting equation for At we have 



At^C. 



C3 
r 



(B13) 



where C2 and C3 are two more integration constants. 
To determine the value of C2 consider the A equation. 



A^ 



T 



-1, 



(B14) 



where Q can be positive or negative. All that is left is to 
determine A^ through the constraint equation (|B14p 



± 1 



2M 




.(B17) 



Phenomenologically i^T < 1 j8|, so the first term in the 
square root is positive. However at large r the second 
term is dominant and is possibly negative. Thus for A^ 
to be real we must satisfy the following bound 



:Q>0. 



(B18) 



To summarise. 
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9 = Diag 



- 1- 



2M Q 



± 1 



2_Q 

K r 

2M Q2 



2M 



^ , 1 + ^ 



,r2,r2 sin2(6») 



K 



(B19) 
(B20) 
(B21) 



Thus the solution is Reissner- Nordstrom (RN) , but where 
the gauge freedom for the field A has been fixed up to a 
sign by the Lagrange multiplier constraint. 

We have obtained a black hole solution in the EF met- 
ric. Its horizons are the ones observed by the gravita- 
tional components g, A, and <j). Any matter fields are 
influenced by the MF metric ^ and it is important to 
consider the MF solution which is the observable frame. 



APPENDIX C: DETAILS OF NON-EXISTENCE 
ARGUMENT FOR BLACK HOLES WITH 
NON-CONSTANT SCALAR 

In this appendix we will argue that there are no black 
hole solutions 'near' to those found in section lTlI Al where 
the scalar is not constant. To support this hypothesis 
we first consider small perturbations to the scalar field 
about this constant scalar black hole solution. We only 
need consider one equation (|C7|) . Since the background 
value for the scalar field is a constant, the equation of 
evolution for the scalar perturbation is simply 



2CiK 



2K{r-2)r+{K-l)(%) -2K^§jr 

(CI) 

Unless Ci = this diverges for two values of r, which 
we denote Tsingi and rsing2 where rsing2 > ''singi- The 
horizon positions are at r+ and r_ with r_|- > It is 
straightforward to show that if i^T > and — 1 < -p^ < 1 
then rsing2 > f+ and so this singularity will occur out- 
side the outermost horizon. A scalar field singularity is 
usually a symptom of a singular geometry. 

Of course, the linear theory breaks down as the scalar 
perturbations become large, and so while it indicates a 
singularity might form outside the horizon if we try to de- 
form the scalar from being constant, it cannot be trusted. 
Hence we also solved numerically the full non-linear the- 
ory with asymptotic data that is close to that of the 
constant scalar RN solution, for particular parameters. 
The full numerical solutions confirm that indeed a naked 
singularity is found, as hinted at by the linear theory, 
with curvature invariants clearly diverging at the singu- 
lar point. We now describe this in detail. 



Using the same metric ansatz (|B3p and the general 
form for the vector field (|B2[) . the r-component of the 
vector field equation is 



Ar[\ + 



87r(0' 



'^2 



kR 



0, 



(C2) 



and so for the case A^. ^Q, this equation determines A to 
be proportional to the square of the proper derivative of 
(f). Substituting A^ from the constraint equation (jB14p . 
and this value of A into the other field equations yields 
Metric (i,<)-component 



2T 

r 



2RT 2TR' lQT:rT{(j)'y 



R 



K 



K 



(C3) 



Metric (r',r)-component 



2 2R 2T' ^ Kr{A'^y 
r r T 

Metric {6 
2rTR' 



167r( 



A'\2 



T 

component 

r^R'T' 



kT 



■ = 0. 
(C4) 



R 



2rT' 



327rr2T(<?!)')2 I&ttt^ Al{cj)'f 



R 



-2Kr^{A'tf 
2r^T" = 0. (C5) 



Vector t-component 



K I -i 

r 

Scalar 



R 



A' 



IGttA 



''^^'^-^2KA'l. (C6) 



r2(2T-A2)' 



(C7) 



where Ci is an integration constant from the second order 
scalar field equation. This is normally associated with a 
scalar mass. Eliminating At and A[ between (jCSp and 
(IC4l) 



167rr((/)' 



/\2 



[TRY 
TR ' 



(C8) 
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and eliminating the same variables between (|C5|) and 
(|C3p . then substituting for cf)' using (|C8p gives 

/2(i?-l) /ri?' \ , r(r')^ 

2T ( A-^ + 1^) + (^ - T ^ T - ' 

(C9) 

which can be integrated once to obtain a first order equa- 
tion in T 

RT{C2 + 4:r^T) = Ar^T^ + Ar^TT' + r^{T'f, (CIO) 



where C2 is a constant of integration. This equation al- 
lows elimination of one of the metric components R from 
the system of equations to numerically integrate. Choos- 
ing to ehminate R and 0' using (jCT)) . from (jUS)) and ((C4| . 
gives the system of equations 



2 /C2T , 9 ^ , ^ ,9 AC^tt(2T + rT'Y\ 
T" = — — + 2 C2 + 2r2T T' - r^T'^ , Cll 



From which (f)' can be obtained during the numerical in- 
tegration through the scalar field equation (jC7p . 

We then integrate inwards from large r, imposing an 
asymptotically flat spacetime, a constant scalar and van- 
ishing Ar component at spatial infinity. Looking first to 
the asymptotic expansion of T and At will allow identifi- 
cation of the free parameters available for the boundary 
condition at large r. Assuming the general asymptotic 
form 



T{r) 
Mr) 



ti 

= to + - + 
r 

r 



U2 



O 



1 



O 



(C13) 
(CM) 



We may use the freedom in rescaling the t coordinate in 
setting to — 1, and the freedom in the r coordinate to set 
ti to —1. Note that this explicitly sets the standard RN 
mass to be positive. It then follows from the constraint 
equation that if we are to have ^ as r — > 00 then 
At —1 and uq = —1. The remaining coefficients may 
be determined by performing a series expansion about 
infinity of the two differential equations (|ClllC12p . We 
find 



T 



At 



1 

h 

r 

-l + ± 



1 



-Go 



4^2 



(C15) 



C2 AnCf 1 



2r 



kK ' 

(C16) 



where the sign choice comes from the sign of the gradient 
of At, equation (|C12p . There is only one consistent choice 
for this, as may be seen by using the Lagrange multiplier 
equation (|B14|) to calculate the corresponding asymptotic 



expansion for A^. To lowest order: 



A^ — 



l-±./-Wl 



-Co 



(C17) 

Thus for the reality of Ar we are forced to choose the 
negative sign. 

We are now left with four parameters: two constants, 
K and K and two free parameters Ci and C2 with an 
additional condition to ensure the reality of At 



l + i^C?_C2>0. 



(C18) 



Additionally, one can identify these parameters with the 
standard Reissner-Nordstrom charge to mass ratio, 



^l-Co 



and a quantity associated with the scalar 
c2 _ ISttC? 



(C19) 



(C20) 



so that the reality condition (|C18[) is now guaranteed. 
Using this definition of 5, the differential equations to 
numerically integrate (|C11IC12[) are reduced to a con- 
stant K, and the two parameters S and Q/M. 

Naively considering the 4> equation of motion (|C7p , one 
sees that a possible divergence occurs when 2T = A^. 
This does occur, and we demonstrate below for the pa- 
rameter choices K — 0.01, k — 0.01, Q/M — 0.01 and a 
small scalar charge, S — 0.001. 

We find that both the EF and MF Ricci scalars diverge, 
and thus it is not a coordinate singularity - this is shown 
in the left panel of figures [HI and [T51 Further, at no point 
outside of the singular position does the gtt component of 
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FIG. 14: EF (solid) and MF (dashed) Ricci Scalars for the 
parameters S = 0.001 and K = n — Q/M = 0.01 as functions 
of the proper distance from the singular surface, illustrating 
a physical singularity. 



In this coordinate system the Einstein equations have 
the usual constraints (the tt and tr components) together 
with a second order evolution equation for i?, the 96 com- 
ponent. We do not use this directly to evolve R, instead 
we use the tr constraint equation itself. Once all vari- 
ables (apart from T) have been successfully evolved to 




0.05 0.1 0.15 

Proper distance from singularity 



the metric vanish - the singularity is not enclosed within 
a horizon, (f), T and At remain finite (though their proper 
gradients diverge) up to the singular point. Interestingly, 
at the singular point grr = R = implying (through 
2T — j4j =0 and = —1) that the radial component of 
the TcVcS vector vanishes at this point. 

We have shown that linear theory suggests no deforma- 
tion of the constant scalar RN solution that has regular 
horizon, and non-trivial scalar. In specific cases we have 
confirmed that the full non-linear theory agrees with the 
linear theory in that attempts to make the scalar be non- 
constant lead to a naked singularity rather than a regu- 
lar horizon. However, we have not shown that there is 
no such solution far in 'solution-space' from the constant 
scalar RN solution. We think it unlikely, although have 
not explored this possibility in detail. 



APPENDIX D: DETAILS OF FULL DYNAMICAL 
NUMERICAL SIMULATION METHOD 

Wo use the following Schwarzschild-like coordinate sys- 
tem for our evolution, 



FIG. 15: Square of the proper derivative of 0, ^ , for the full 

integration (solid) and for the linearised theory (dashed). A 
larger S would see a larger discrepancy between the singular 
positions for the full and the linearised theory. This plot also 
demonstrates the agreement of the two curves away from the 
singularity. 



-T^{t, r)df + e^^^'^Ur^ + r^d^l (Dl) 



the next spatial slice, we integrate across the grid in the 
radial direction to obtain T on the slice using the rr equa- 
tion. A simplification is made for this radial integration 
of T; we use the value of T from the old spatial slice to 
compute the contribution to the right hand side of the 
rr-Einstein equation, for the sake of computational run 
time. Using this approximation, the rr component is a 
first order differential equation in T. 

We find that the most stable way to evolve the vector 
field is to evolve A, rather than evolving A and then calcu- 
lating A through contraction of the vector field equation. 
The evolution equation we use for A is given by taking 
the divergence of the vector equation 



VcV^F'^" + V„(AA") + Va{<J^A^cl>,0 g'^'^ct^,^ ) 



(1 - e-^'^)5"^A^i; 



(D2) 
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The first term vanishes through antisymmetry of F and 
the symmetry of the Ricci Tensor. The third term will 
contain second time derivatives of (j), and so we substitute 
for this using the </> equation of motion. The last term 
will contain second derivatives of the x field. For this 
we do not substitute from the x field equation as we are 
able to use the conservation properties of T to obtain 
a simpler (and so numerically advantageous) expression. 
Consider the definition of T from a variation of the action 



Specifically, consider the case of a diffeomorphism gener- 
ated by the vector field itself, 

6g'^'^ = £a9°''^ = V^"^''), (D5) 
6A^' = fAA" ^ 0, (D6) 
6(1) = £a<I) = A"V„^, (D7) 

so that 

6 A, = 6{g,^A'') = -A^g^pSg'^P = -A„g,^V("A«. 

(D8) 

Inserting all of this into (jD3|l and integrating by parts 
from the terms containing variations of the inverse EF 
metric, we obtain the following expression 

A^V^fo^f} = + (1 + e^'*'*) A'^AP) A^'V ^.cf), (D9) 

which we can then use instead of the x evolution equation 
in the evolution of A. 

Thus both the metric components and A are evolved 
with first order differential equations, whilst A, and x 
are evolved at second order. The origin boundary condi- 
tions are <j),r = 0, A'' = 0, T = 1, i? = 0,r,^ = 0, R,r = 
and Xir — 0. 

For the scalar shell collapse in section (jIV A[) the initial 
conditions on the t — Q spatial slice, are (j) — (pc, A^ — 0, 
R,t — T,t ~ 4>,t = Xit ~ 0. We consider two simulations 
with differing initial conditions for -F)t,y, and these are 
discussed in the text. The initial R and T configurations 
are specified by solving the rr and tt Einstein equations 
for a choice of Gaussian initial data on x- The tr-Einstein 
equation is automatically satisfied for this initial data. 

For the boson star in section pVBp the initial con- 
ditions for all fields are fixed by requiring for a static 
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with respect to the MF metric, 

SS = -^f^0^6r^, (D3) 

and one may rewrite the variation of g in terms of vari- 
ations of the other fields, through the disformal relation 
Misee i) 



(D4) 

I 

star configuration initially. This configuration was found 
using a radial shooting method, as discussed in the text. 

Second order finite differencing was used to discretize 
the equations of motion. For the scalar-shell collapse 
simulation, we used resolutions up to Ar = 0.05 and 
At — 0.0005, whilst for the Boson star simulation we 
used resolutions up to Ar = 0.05 and At = 0.00005. 
Each simulation had 400 simulation sites per slice, and 
took on the order of a day to run using an average desktop 
computer. 



APPENDIX E: CONSTRAINT TESTING 

We ran simulations at a variety of resolutions to check 
convergence which was seen in accord with our second 
order finite differencing. To test the constraint equations 
(and convergence), we consider the black hole dynamics 
for some representative initial data as in section IIV Al 
and take the absolute value of the constraint equation 
(l.h.s. - r.h.s.) at a grid point, and then average this over 
all grid points (labelled by i) in a large physical area. We 
then compute this average for various spatial and tem- 
poral resolutions, keeping the physical region fixed. We 
denote the set of all grid points in this region as a. Figure 
1161 shows this sum for the tt component of the Einstein 
equation against logj^g fo'' four spatial steps Ar — 0.2, 
0.1, 0.05 and 0.025. We clearly see agreement with ex- 
pected constraint behaviour for second order differencing 
as the spatial and temporal resolutions are reduced. Sim- 
ilar checks were performed for the boson star simulation 
of section HVB 
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FIG. 16: The average of the tt constraint equation over simu- 
lated grid points for the black hole of section IIV Al for differ- 
ent time steps. Each curve corresponds to a different spatial 
resolution. 
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